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u  abstract 


This  work,  presents  the  design  considerations  and  the  calibration  of  a  shock 
tunnel  for  aerodynamic  measurements  at  hypersonic  Mach  numbers.  The  shock  tunnel 
is  composed  of  a  shock  tube  which  is  used  for  the  production  of  high  temperature 
and  high  pressure  air  reservoir,  which  is  then  expanded  through  a  supersonic  nozzle, 
thereby  attaining  high  Mach  number  flow  in  the  test  section. 

In  order  to  obtain  maximum  test  time  .the  "tailored”  mode  of  operation  is  used 
for  this  shock  tunnel.  A  method  for  calculating  the  optimum  ratio  of  the  compression 
chamber  length  to  the  low  pressure  tube  length  for  a  given  total  shock  tube  length 
that  will  result  in  a  maximum  test  duration  is  presented.  Shock  tube  experimental 
results  at  "tailored"  condition  for  hydrogen  (driver  gas)  and  air  (as  the  Jriven 
gas),  both  at  room  temperature,  indicate  that  by  proper  design  a  reservoir  of  air 
with  an  enthalpy  of  2000  Btu/lb.  can  be  produced  for  a  duration  of  mor.e  than  2.0 
msec,  for  a  total  shock  tube  length  of  9  meters. 


The  flow  in  the  test  section  is  calibrated  using  a  total  pressure  rake  and  by 
static  pressure  measurements  and  measurements  of  stagnation  point  heat  transfer 
rates  at  various  points  in  the  test  section.  The  flow  conditions  are  then  calculated 
using  the  theoretical  relations  for  stagnation  point  heat  transfer  rates.  Thg  shock 
tunnel  test  section  has  been  thus  calibrated  at  two  nozzle  positions:  one  A/A  ■  250 
resulting  in  a  nominal  Mach  number  of  5.5  ami  A/A*  «  500  resulting  in  a  nominal  Mach 
number  of  7. 


The  experimental  results  show  a  reasonably  uniform  flow  in  the  horizontal  plane 
with  a  relatively  thin  boundary  layer.  While  in  the  vertical  plane  only  a  central 
core  of  about  half  test  section  height  has  uniform  flow  conditions  and  is  suitable 
for  aerodynamic  testing. 
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I .  INTRODUCTION 


The  shock  tunnel  utilizes  for  its  operation  the  reservoir  of  high 
temperature,  high  pressure  air  at  the  end  of  the  shock  tube.  This  air  is 
then  expanded  through  the  hypersonic  nozzle  to  high  Mach  numbers.  The  flow 
regions  in  the  shock  tube  are  shown  in  a  schematic  diagram  in  Figure  1.  The 
hot  gas  enclosed  between  t.he  first  shock  wave  and  the  contact  surface  (region 
2)  may  be  utilized  for  aerodynamic  measurements.  The  duration  of  this  flow 
depends  upon  the  shock  Mach  number  and  the  length  of  the  low  pressure  section. 

The  time  duration  for  a  unit  length  of  low  pressure  section  as  a  function  of  • 
the  shock  wave  Mach  number  is  presented  in  Figure  2.  If  now,  a  supersonic 
nozzle  is  connected  to  the  end  of  the  shock  tube  it  is  possible  to  accelerate 
the  shock  heated  air  in  region  2  to  a  hypersonic  Mach  number.  This  arrange¬ 
ment  of  the  "flow  ihrough"  shock  tunnel  is  limited  by  the  very  short  test 
duration.  As  shown  in  Fig.  2,  this  time  can  be  few  hundreds  of  microseconds  for 

a  tunnel  of  acceptable  length,  furthermore,  this  test  duration  decreases  with 
increasing  initial  shock  Mach  number.  Another  method  of  operating  the  shock 
tunnel  utilizes  the  hot  gas  behind  the  reflected  shock  wave  at  the  end  of  the 
shock  tube  at  the  conditions  of  "tailored  interface".  This  condition  is  achieved 
when  the  shock  wave  which  is  reflected  from  the  end  of  the  shock  tube  encounters 
the  contact  surface  without  any  reflected  disturbances.  At  "tailored"  conditions 
a  much  longer  (few  milliseconds)  test  time  is  available  for  shock  tunnel  applications. 
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II.  CALCULATION  OF  THE  TAILORED  SHOCK  TUBE  CHARACTERISTICS 


2.1.  Calculation  of  Tailoring  Conditions 


The  calculation  of  the  tailored  conditions  is  described  in  Ref.  3.  A 
short  review  of  this  calculation  is  presented  in  order  to  facilitate  the 
description  of  the  present  study 

The  conditions  which  must  exist  at  the  point  of  encounter  between  the 
reflected  shock  wave  and  the  contact  surface  is  that  the  static  pressure  and 
velocities  must  be  equal  on  both  sides  of  the  contact  surface  i.e.  (see  Fig. 


U6*U7 

P6*P7 

The  equation  which  relates  the  velocities  and  pressures  on  the  two  sides 


of  a  moving  shock  wave  is 


where: 


-  upstream  of  the  shock  wave 

-  downstream  of  the  shock  wave 


t  _  r  2  1 0 .5 

{  •  'tftt1 


y  ♦  1 

Y  -  1 


and  the  positive  velocity  direction  is  in  the  direction  of  movement  of  the 


shock  wave.  From  these: 
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Since  U5  *  0 


{4<P73 


»l°4P73  *  »"°'5 


(2) 


VP65  -  ‘>(“lP6S  *  « 


(3) 


For  the  case  in  which  the  reflected  disturbances  are  rarefaction  waves  wc 
obtain: 


where: 


is  1 

±  =  -  — —  (i  - 


(4) 


Since  »  0  is  the  contact  surface  velocity  after  the  encounter,  it  follows 
from  Eqs .  3  and  4  that  the  contact  surface  moves  toward  the  nozzle  when  the 
reflected  disturbances  are  shock  waves  and  away  from  the  nozzle  when  the  re¬ 
flected  disturbances  are  rarefaction  waves. 


From  the  conditions  *  U^,  we  obtain  from  Eqs.  2  and  3  for  the  re¬ 
flected  shock  wave: 
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and  for  the  reflected  wave  (from  Equations  2  and  4) 
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From  the  condition  P, 
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7,  we  obtain 
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When  the  reflected  diturbances  from  the  contact  surface  is  a  Mach  wave, that 
is  to  say  P&5  =  1,  we  obtain  from  Eq.  (7)  that  P^^  «*  P^.  In  this  case 
from  Eqs.  (5)  and  (6)  (since  P2  8  P3  )  one  obtains 
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but  from  Eq.  (1)  we  may  write 
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when  Uj.  «  0  (U2  is  in  the  direction  opposite  to  the  direction  of  the  shock 

wavej,  and,  from  Eqs.  (8),  (9)  end  since  U2  «  we  obtain: 
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Eq.  10  presents  the  conditions  for  not  producing  reflected  disturbances  from 
the  contact  surface.  This  equation  gives  the  relation  between  a^  and  Pg2« 

As  P52  is  a  function  of  P41,  Eq.  (10)  can  be  used  to  obtain  the  relation 

between  a0_  and  P..  for  "tailored1*  condition.  In  the  shock  tube  flow  one 
*3  4 1 

can  derive  the  following  relation: 
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Since  ?21  is  a  function  of  P41,  an  additional  equation  is  obtained  which 
gives  a  relation  between  a23  and  P^. 

Thus  Eqs.  (10)  and  (11)  define  the  conditions  for  no  disturbances  du  •  to 
interaction  of  the  reflected  shock  wave  and  the  contact  surface.  Equation  (11) 
presents  the  conditions  required  in  the  shock  tube  in  order  to  obtain  the  "tailored" 
value  for  a23,  for  a  given  combination  of  driver  and  driven  gases  at  given  initial 
temperatures  (the  temperature  influences  the  quantity  a41).  The  variation  of  a23 
with  Ms  for  air  (  as  the  driven  gas)  and  various  driver  gases  is  presented  in 
Fig.  3.  Thus,  at  this  "tailored"  conditions,  we  obtain  between  the  contact  surface 
and  the  end  of  the  shock  tube  a  reservoir  of  high  temperature,  high  pressure  air 
of  uniform  properties  and  free  of  disturbances.  By  allowing  this  air  to  flow  through 
a  convergent-divergent  nozzle  we  obtain  high  Mach  number  flow  in  the  test  section 
of  the  shock  tunnel.  This  is  done  by  attaching  a  nozzle  of  relatively  small  throat 
area  (in  this  case,  4%  of  the  cross  section  area),  at  the  end  of  the  shock  tube. 

In  this  case  a  complete  reflection  of  the  incident  shock  wave  still  occurs  before 
the  shock  heated  air  is  expanded  through  the  nozzle. 

It  should  be  pointed  out  that  at  conditions  slightly  different  from  the 
tailoring  conditions,  the  reflected  disturbances  from  the  contact  surface  are 
weak  and  there  still  exists  a  reservoir  of  heated  air  whose  conditions  are  very 
nearly  uniform.  It  is  pointed  out  in  Ref. [4]  that  for  5.7  <  Mg  <  6.3,  the 
deviation  from  the  tailoring  pressure  is  at  most  *  10%.  A  theoretical  develp- 
ment  of  disturbance  strength  for  the  case  of  deviations  from  the  tailoring  con¬ 
ditions  may  be  found  in  Ref.  [5]. 
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2.2.  Calculation  of  the  Duration  of  Flow  at  Tailored  Conditions 

The  duration  of  the  experiment  in  the  test  section  of  the  nozzle  is 
limited  by  the  following  factors:  The  depletion  of  the  reservoir  of  the  hot 
gas  (i.e.  the  arrival  of  the  contact  surface  to  the  entrance  of  the  nozzle), 
and  due  to  the  arrival  of  the  rarefaction  wave  to  the  nozzle  entrance. 

a)  Calculation  of  the  Mass  Depletion  Time 

The  mass  depletion  time  can  be  calculated  equating  the  mass  flowing  out 
through  the  nozzle  to  the  mass  depletion  due  to  the  velocity  with  which 
the  contact  surface  is  carried  towards  the  nozzle  entrance,  that  is 

>tVt  *  <Wc.  (12) 

where 

U  is  the  contact  surface  velocity  towards  the  nozzle. 

Since  Uj/aT  «  1,  then 


p/t  ■ 

:  r- PcAcu„ 
a^  5  5  c .  s . 

It  follows  that: 

Uc.s. 

Aypy 

87  Vs 

(13) 

aT 

The  ratios  —  and 
a5 

PT 

—  can  be  calculated  by  assuming  that 
P5 

the  flow 

up  to  the  throat  is 

isentropic.  Then, 

Uc,s. 

1  +  y 

.  ,  .1 . 

a5  Ac  1  2  } 

(14) 

For  the  case  in  which  the  air  between  the  contact  surface  and  the  end 
of  the  shock  tube  deviates  from  an  ideal  gas  then  the  ratio  as/Lp  and 
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and  Pg/'j  way  be  obtained  from  charts  for  air  at  thermodynamic  equilibrium. 


rhe  depletion  time,  can  be  now  estimated  to  be, 

rwh 


c.s. 
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where, 


t  -  is  the  distance  of  the  contact  surface  from  the  end 
c.s. 

of  the  shock  tube  and  is  calculated  from  the  relation 
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In  this  case 
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For  specific  initial  conditions,  for  which  we  wish  to  obtain  "tailoring'* 
(i,e.  known  y,  A_,  a,.,  U,,  U  and  U  ),  t  _  will  vary  with  the  ratio 
Ag/AT  and  with  L  as  indicated  in  Eq.  (18).  Of  course  this  time  is  limited  by 
the  maximum  length  of  low  pressure  section  which  is  limited  by  boundary  layer 
growth . 


b)  Calculation  of  the  Time  for  the  Arrival  of  the  Rarefaction  Wave  to  the 
Rozzlo", 


The  rarefaction  wave  travels  with  the  local  speed  of  sound  and  its  propagation 
can  be  calculated  by  the  procedures  described  in  Ref.  1. 


2.3.  Calculation  of  Optimum  llTailoredllShock  Tube  Design 


We  will  define  an  optimum  design  for  a  shock  tube  that  will  produce  the 
longest  test  flow  duration  in  a  shock  tunnel  at  "tailored"  conditions  when  the 
total  length  of  the  shock  tube  is  given.  This  optimum  design  will  exist  when 
the  shock  tube  is  so  proportioned  so  that  the  wave  head,  the  rarefaction  wave 
tail  and  the  contact  surface  arrive  simultaneously  at  the  nozzle  entrance 
position.  This  condition  is  presenred  in  the  following  relation: 


+  L*  ♦  L 
a4  a  +  0 


+  B  x  C  x  L  (19) 


time  in 
which 
the 
rare¬ 
faction 
wave 
tail 
arrives 


the  time  in 
which  the 
first  shock 
wave  arrives 
at  the  nozzle 

A 


Lost  duration 
due  to  the 
depletion  of 
the  reservoir 

A 


time  in  which 
rarefaction 
wave  head 
arrives 


the  time 
in  which 
the  reflected 
shock  wave 
travels  the 
distance 
between  nozzle 
and  the  point 
of  encounter 
with  the  contact 
surface 


the  time  in  which  the  contact  surface  arrives. 


where: 


0  -  the  average  velocity  at  which  the  reflected  rarefaction  wave 
passes  through  the  expansion  fan. 
a  -  average  sonic  velocity  in  the  expansion  fan 
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Thus,  in  order  to  obtain  an  optimum  shock  tube  we  require  that  the  follow¬ 
ing  relation  holds  for  the  optimum  L, 
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and  in  this  case  the  appropriate  value  of  Ag/A^,  is: 
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The  additional  relation  which  follows  from  Eq.  (19)  is 
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-1).  (2D 


(22) 


Substituting  the  value  of  Ag/Aj,  obtained  from  Eq.  (21),  in  Eq.  (22)  we 
obtain  the  relation  for  optimal  lengths  of  L  and  l'  "  be. 


DL'  *  (rr^-s - L’  (23) 

us  ‘  a3  a  +  0 

The  application  of  the  present  calculation  to  the  shock  tube  used  in 
this  investigation  indicates  the  proper  lengths  that  should  be  used  in  order 
to  obtain  optimum  test  flow  duration  and  are  shown  in  Figs.  4  and  S. 
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2.4.  The  Additional  Reservoir  Pressure 


In  addition  to  the  wave  system  produced  with  the  bursting  of  the  first 
diaphragm  which  was  described  previously,  there  is  a  supplementary  shock  wave 
that  reaches  the  region  of  standing  air  (region  5)  and  causes  a  pressure  rise 
in  this  region.  It  should  be  pointed  out  that  this  shock  wave  appears  also 
when  "tailoring"  conditions  exist.  The  reason  for  the  appearance  of  this  shock 
wave  was  discussed,  based  on  theoretical  considerations,  in  reference  (6).  This 
theory  will  be  briefly  reviewed  herein.  In  the  shock  tube  driver  system  present¬ 
ed  in  Fig.  6  the  bursting  of  the  diaphragm  due  to  the  high  pressure  in  region  4 
will  produce  an  expansion  wave  which  will  advance  into  region  4.  This  expansion 
wave  causes  an  unsteady  expansion  until  a  new  condition  (designated  as  3u)  is 
attained.  This  condition  will  depend  on  the  initial  conditions  of  the  driver 
and  driven  gases.  The  equation  relating  the  pressure  and  velocity  in  this  case 

will  be  (relation  between  the  quantities  in  an  expansion  wave): 

2Y„ 


?3u  „  Y4  "  1  U34,y4  ’  1 

P7  *  (1-—  7? 


(24) 


The  expansion  of  the  driver  gas  through  the  nozzle  at  the  diaphragm  section  into 

region  1  is  calculated  by  the  use  of  the  steady  flow  equations  through  the  nozzle. 

The  relation  between  the  pressure  and  the  velocity  is  then: 

y4 


r'-c  Y4  ■  1  U35  2  Vl 

pf"  »  -  'sf1  > 
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(25) 


where  is  the  condition  of  the  gas  following  the  steady  expansion.  From  Eqs. 

(24)  and  (25)  we  see  that  if  we  start  from  the  same  initial  conditions  and  attain 
the  same  velocity,  the  resulting  pressures  P3u  and  P3s  will  be  different.  There 
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is  only  one  velocity  at  which  these  pressures  are  equal.  This  velocity  can 
be  calculated  by  equating  the  pressures  in  Eqs .  (24)  and  (25) . 

This  condition  is  called  in  Ref.  (6)  the  "cross  over  point".  Since  this 
condition  does  not  correspond  to  the  "tailored"  condition  it  is  required  to 
introduce  an  additional  unsteady  wave  which  provides  the  matching  boundary  con¬ 
ditions  in  the  "tailored"  shock  tube.  This  patching  process  may  be  accomplished 
by  either  a  shock  wave  or  a  rarefaction  wave  moving  to  the  right  (in  the  direction 
of  the  supersonic  nozzle)  or  to  the  left  respectively. 

The  conditions  which  govern  the  appearance  of  these  waves  and  their  direction 
depend  upon  two  factors.  One  condition  is  the  ratio  between  U_  or  M  for 
tailoring  and  for  a  cross  over  point  and  the  second  is  the  area  ratio  of  the  nozzle. 
Reference  (6)  gives  details  for  various  cases.  Those  of  interest  for  this  work  are 
the  following:  The  case  where  M„/M  <  1,  as  is  illustrated  in  the  _  U  dia- 
gram  of  Fig.  7.  Presented  in  this  figure  are  the  curves  for  steady  and  unsteady 
expansions.  The  point  A  on  the  unsteady  expansion  curve  corresponds  to  the 
velocity  and  pressure  obtained  at  the  tailoring  conditions.  The  corresponding 
point  on  the  steady  expansion  line  depends  on  P^  and  on  the  area  ratio  of  the 
nozzle.  In  case  where  there  is  no  nozzle  at  the  diaphragm  section,  that  is  6*1, 
patching  is  achieved  by  the  expansion  wave  moving  to  the  left  and  shock  wave  moving 
io  the  right.  This  shock  wave  reaches  the  resevoir  region  (region  5)  and  raises  the 
pressure  in  this  region.  This  phenomenon  is  well  illustrated  in  the  experimental 
results  obtained  in  this  work.  When  6  becomes  less  t..un  1,  the  strength  of  the 
two  waves  decreases,  and  when  6  *  >e  patched  conditions  are  obtained  by  only  one 
shock  wave  which  moves  to  the  righ  ,  An  additional  decrease  in  6  results  in  two 
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shock  waves,  one  of  which  will  move  to  the  right  and  the  second  to  the  left. 
However,  the  shock  wave  which  moves  to  the  right  will  become  weaker  with  de¬ 
creasing  6  ,  while  the  one  which  moves  to  the  left  with  increase  in  strength. 

When  6  =  ^jeai,  patching  conditions  will  be  provided  by  only  one  shock  wave 
which  moves  to  the  left.  It  is  clear  that  in  this  case  the  air  behind  the 
reflected  shock  at  the  right  end  of  the  shock  tube  will  not  be  disturbed. 

Additional  reduction  of  6  causes  a  leftwards  moving  shock  wave  and  a  right- 
ward  moving  rarefaction  wave.  This  rarefaction  wave  will  reduce  the  pressure 
of  the  air  in  regi  t  S. 

Therefore,  it  is  clear  that  if  we  wish  to  prevent  the  additional  shock  wave, 
we  must  install  an  additional  nozzle  in  region  4.  A  method  for  the  calculation  of 
^ideal  *s  Presented  in  Reference  (6).  In  the  present  shock  tube  using  H2/Air 
combination  and  »  6.0  the  required  value  for  6  is  0.382.  This  condition  was 
well  verified  in  the  experimental  program  and  is  discussed  .'u  Part  III  (  1  ). 

The  installation  of  a  nozzle  in  region  4  removes  the  head  and  the  tail  of  the 

rarefaction  wave  which  is  produced  at  the  bursting  of  the  diaphragm.  However,  the 
collision  between  the  first  reflected  shock  wave  and  the  leftward  moving  patching 
wave  produces  a  reflected  disturbance  in  the  form  of  a  rarefaction  wave  which  moves 
to  the  right.  The  arrival  of  this  rarefaction  wave  at  region  5  will  terminate  the 
useful  test  duration.  It  is  worthwhile  to  note  that  the  experimental  results  of 
this  work  show  that  this  rarefaction  wave  arrive  at  region  5  at  about  the  same  time 

that  the  first  rarefaction  wave  head  would  arrive  if  there  were  no  nozzle.  It  is 

important  to  consider  the  effect  of  this  additional  w*.ve  on  the  shock  tube  testing. 
There  are  cases  where  it  was  suggested  to  incorporate  a  nozzle  at  the  diaphragm  in 
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order  to  stop  the  rarefaction  wave  (Ref.  7),  and  therefore  eliminate  the 
propagation  of  this  patched  wave  into  region  5.  However  this  additional 
wave  can  also  be  used  to  increase  the  pressure  in  region  (5)  for  certain 
shock  tube  geometry  (as  was  done  in  the  present  investigation)  at  a  certain 
loss  in  test  duration. 
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III.  FLOW  IN  THE  HYPERSONIC  SHOCK  TUNNEL  NOZZLE 


3.1.  Calculation  of  the  Stagnation  Conditions  (Conditions  in  Region  5) 

From  knowledge  of  and  conditions  in  region  1  of  the  shock  tube 
it  is  possible  to  calculate  the  conditions  behind  the  reflected  shock  wave 
(region  S)  which  constitute  the  stagnant  air  reservoir  for  the  flow  in  the 
supersonic  nozzle.  Agreement  between  the  calculated  quantities  and  the 
actual  values  was  verified  from  recordings  of  the  pressure  at  the  end  of  the 
shock  tube;  the  experimental  results  deviate  only  slightly  from  those  obtained 
in  the  theoretical  calculations.  The  thermodynamic  properties  (assuming 
equilibrium  conditions)  were  obtained  from  References  8  and  9.  These  calculated 
quantities  are  correct  for  the  case  in  which  the  patching  wave  is  removed  by 
use  of  a  nozzle  in  region  4  (discussed  in  Section  2.4).  In  the  case  in  which 
the  patching  wave  is  not  removed,  it  was  found  experimentally  that  the  final 
pressure  region  5  was  larger  by  factor  of  1.6  than  that  in  a  system  in  which 
the  patching  wave  is  removed.  It  should  be  remembered  that  this  pressure 
rise  results  from  two  waves,  the  patching  wave  which  passes  through  region 
5  and  the  reflection  of  this  wave  from  end  of  the  shock  tube.  If  we  assume 
that  all  of  the  rise  in  pressure  results  from  one  wave  only  then  for  this 
pressure  rise  the  corresponding  rise  in  the  total  enthalpy  and  in  the 
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stagnation  temperature  is  only  10%.  Therefore,  an  average  value  for  the 
total  enthalpy  may  be  taken  in  the  reservoir  with  or  without  patching 
wave  with  reasonable  accuracy.  Of  course  the  pressure  and  the  density 
will  be  different  in  the  two  cases.  Table  i  presents  typical  shock  tube 
conditions  for  the  operation  of  the  tailored  shock  tunnel.  At  these  con¬ 
ditions  there  exists  (Fig.  8,  which  was  taken  from  Reference  (10))only  one 
significant  reaction  in  the  air,  i.e.  the  dissociation  of  the  oxygen  molecules 
into  atoms. 

The  steady  flow  in  the  supersonic  nozzle  (after  the  flow  starting  process) 
can  be  assumed  to  be  isentropic,  while  the  flow  may  be  either  in  thermodynamic 
equilibrium  or  fronzen  conditions. 

A  computer  chart  for  the  calculation  of  the  equilibrium  flow  knowing  the 
reservoir  conditions  is  presented  in  Fig.  9.  Results  for  isentropic  thermo¬ 
dynamic  equilibrium  flow  in  hypersonic  nozzle  are  presented  in  Ref.  11. 
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3.2.  The  Starting  Process  in  the  Shock  Tunnel  Nozzle 

Prior  to  the  start  of  the  experiment  a  cellophane  diaphragm  is  inserted 
between  the  shock  tube  and  the  hypersonic  nozzle.  The  air  pressure  in  the 
nozzle  is  lower  than  that  in  the  shock  tube.  It  is  clear  that  if  the  initial 
pressure  in  the  nozzle  will  be  higher  than  the  static  pressure  of  the  steady 
hypersonic  flow  then  the  proper  flow  in  the  nozzle  cannot  be  established 
(Ref.  15). 

The  flow  immediately  after  bursting  of  the  diaphragm  is  similar  to  the 
flow  immediately  after  bursting  of  the  diaphragm  in  the  shock  tube.  Here  also, 
due  to  the  pressure  difference  across  diaphragm,  a  shock  wave  is  produced  which 
travels  to  the  right  (in  the  direction  of  the  dump  tank)  and  a  rarefaction  wave 
which  moves  to  the  left  (in  the  direction  of  the  shock  tube) .  This  wave  structure 
associated  with  the  flow  starting  phenomena  in  the  hypersonic  nozzle  is  very  com¬ 
plicated  due  to  the  variation  in  the  cross-sectional  area  of  the  nozzle  and  due  to 
the  secondary  disturbances  produced  following  interaction  between  the  shock  wave 
and  the  rarefaction  wave.  Fortunately,  quantitative  evaluation  is  not  of  interest 
since  the  main  requirement  is  that  these  disturbances  will  be  swept  downstream  as 
quickly  as  possible.  Experimental  results  in  Refs.  15  and  16  show  that  the  greater 
the  initial  pressure  ratio  between  the  shock  tube  and  the  hypersonic  nozzle  the 
shorter  is  the  time  during  which  these  disturbances  are  swept  away. 

Flow  starting  disturbances  which  occur  in  the  shock  tube  may  be  qualitative¬ 
ly  explained  with  the  aid  of  a  p  -  U  diagram  similar  to  that  used  for  the  explan¬ 
ation  of  the  patching  wave  disturbances  in  the  shock  tube  (Section  2.4).  It  is  not 
possible  to  connect  the  initial  pressure  in  the  nozzle  and  the  static  pressure  of 
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the  steady  flow  by  a  rightward  moving  shock  wave  and  a  leftward  moving  rare¬ 
faction  wave  only.  An  additional  unsteady  phenomenon  is  required,  a  rare¬ 
faction  wave  or  shock  wave.  When  the  initial  pressure  in  the  shock  tunnel 
nozzle  is  low,  patching  between  the  steady  and  the  unsteady  flows  is  accommplish- 
ed  by  a  leftward  moving  rarefaction  wave  and  a  rightward  moving  shock  wave.  As 
the  initial  pressure  in  the  nozzle  is  increased  the  leftward  moving  disturbance 
is  transformed  from  the  rarefaction  wave  to  a  shock  wave.  It  should  be  noted  that 
while  these  disturbances  move  with  respect  to  each  other,  they  are  all  swept  right¬ 
wards  by  the  air  flowing  in  this  direction  at  hypersonic  speed.  The  weaker  these 
disturbances  are  and  the  higher  the  flow  velocity  is,  these  disturbances  are  then 
swept  faster  passed  the  test  section.  The  duration  of  the  passage  of  these 
disturbances  of  the  nozzle  beyond  test  section  is  actually  the  lost  test  time 
due  to  the  starting  process. 

In  the  case  of  frozen  flow,  calculations  for  the  flow  parameters  at  various 
reservoir  conditions  and  various  Mach  numbers  for  which  the  flow  freezes  are 
presented  in  Ref.  12. 


3.3.  The  Calculation  of  the  Stagnation  Point  Heat  Transfer  in  the  Hypersonic 
Shock  Tunnel  Flow. 


The  stagnation  point  heat  transfer  in  hypersonic  flow  has  been  calculated 
in  References  13  and  14,  and  can  be  presented  in  the  following  equation: 
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where: 


For  frozen  flow 

For  equilibrium  flow 

and  h„  *  h?a 
c  As 


-  c  *  0.63. 

-  c  =  0.52. 


Since  L  is  close  to  1,  the  terms  h  /h  is  less  than  one  and  o  is 

c  s  s 

small  due  to  the  relatively  high  pressure  at  the  stagnation  point,  the 
equation  may  be  written  in  the  form: 


-0.6. 


0.5  ,VV 0.1  CS.(h„  -  h) 


dU, 


’w 


0.765  Pr”'"v'(p  y  y"'  ( -^-)u“  ^sV"S  w' 


(27) 


S"S' 


P  u 
HSMS 


The  velocity  gradient  at  the  stagnation  point  of  a  sphere  will  be: 
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The  value  of  A  is  1.41  for  Newtonian  pressure  distribution  on  the  sphere  and 
p  is  then  the  static  pressure  of  the  undisturbed  flow.  If  the  flow  about  the 
sphere  is  assumed  to  be  like  the  incompressible  potential  solution  then  A  *  2.12 
and  p  is  the  static  pressure  behind  the  detached  shock  wave.  The  evaluation 
of  these  values  is  presented  in  Appendix  A. 

For  small  values  of  h  /h  ,  which  is  the  situation  in  shock  tunnel  (since 

w  s 

the  model  is  at  room  temperature)  then 
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n  ,  n  r  p  0.25  p  0,25 

Let  K  *  0.763  Pr  ,0  Au  ,  and  since  (-i)  =  (-1)  ,  then 
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the  following  equation  is  obtained  (based  on  the  known  conditions  in  the  shock 
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The  quantity  qw  -g-y  which  is  calculated  from  Eq.  (31)  for  air  in  thermo- 
P1 

dynamic  equilibrium  as  a  function  of  the  flow  Mach  number  in  the  shock  tunnel 


at  various  stagnation  enthalpy  levels  is  shown  in  Fig.  10. 
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IV.  THE  experimental  system 

4.1.  The  Shock  Tunnel 

The  main  part  of  the  shock  tunnel  and  its  overall  dimensions  are  shown 
in  Fig.  11.  The  shock  tunnel  consists  of  two  main  parts:  the  shock  tube 
and  the  hypersonic  nozzle.  Photographs  of  the  shock  tunnel  are  presented  in 
Fig.  12.  The  high  pressure  section  of  the  shock  tube  consists  of  a  circular 
cross-section  with  an  internal  diameter  of  80mm.  and  an  external  diameter  of 
115  mm.  and  is  designed  for  operation  to  internal  pressures  of  200  atm.  This 
tube  is  made  of  steel  and  is  2  meters  long.  One  end  of  the  tube  is  closed  by 
a  flange  to  which  the  driver  gas  feed  system  is  connected.  The  other  end  is 
connected  to  the  low  pressure  channel  and  at  this  end  the  diaphragm  is  also 
located.  The  low  pressure  channel  is  7  meters  long  with  a  square  75  x  75  mm. 
cross-section.  Although  this  low  pressure  section  is  limited  in  the  pressure 
range  of  its  operation,  the  two  dimensional  flow  present  in  the  square  cross- 
section  is  preferred  for  aerodynamic  studies  when  it  is  used  as  a  shock  tube. 

A  test  section  is  located  near  the  end  of  this  tube  with  windows  through  which 
it  is  possible  to  install  models  for  research  in  the  shock  tube  flow.  Various 
measurements  can  be  made  with  transducers  installed  through  several  ports  which 
are  placed  along  the  last  three  meters  of  the  shock  tube.  A  convergent  divergent 
hypersonic  nozzle  is  installed  at  the  end  of  the  shock  tube.  The  convergent 
section  reduces  the  cross-section  area  from  the  75  x  75mm.  channel  to  a  throat 
of  15  x  15  mm.  Downstream  of  the  throat  provisions  are  made  for  holding  the 
cellophane  diaphragms  which  separate  the  shock  tube  initial  pressure  from  that 
of  the  hypersonic  nozzle.  The  divergent  section  is  composed  of  two  sections. 


22 


* 

In  the  first  section,  the  nozzle  diverges  to  the  area  ratio,  A/A  *  78.  This 
section  starts  in  a  horizontal  position  and  ends  with  a  deflection  of  10  degrees 
through  a  deflection  plate.This  deflector  plate  is  introduced  in  order  to  centi- 
fuge  out  the  cellophane  particle:  which  are  carried  with  the  flow  and  thereby 
reduce  the  damage  to  models  in  the  test  section.  It  is  worth  noting  that  this 
arrangement  was  found  to  be  efficient  since  heat  transfer  gauges,  which  are 
very  sensitive  to  particle  damage,  were  almost  undamaged  during  the  progress 
of  these  experiments.  The  second  section  of  the  nozzle,  which  protrudes  slight¬ 
ly  into  the  first  .section,  starts  with  a  cross-sectional  area  smaller  than 

that  at  the  end  of  the  first  section;  this  permits  the  removal  of  the  boundary 

* 

layer  of  the  first  section  and  makes  it  possible  to  obtain  an  area  ratio,  A/A 
larger  than  would  be  possible  in  a  one  Stage  nozzle  .  The  width  of  this 
second  section  is  constant  at  200  mm.  while  the  height  is  varied  by  a  wedge  type 
nozzle.  The  wedge  angle  can  be  varied,  so  that  it  is  possible  to  vary  the 
entrance  area  to  the  test  section  and  therefore  the  flow  properties  in  this 
section.  Fig.  13  shows  the  second  section  of  the  nozzle  (with  the  first  section 
removed)  and  the  associated  equipment  including  the  0.25  m  dump  tank  which  is 
connected  to  the  end  of  the  hypersonic  nozzle,  the  boundary  layer  removal  system, 
the  diffusion  pump  etc.  The  test  section  is  incorporated  into  the  second 
section  of  the  nozzle,  and  a  model  support  strut  is  installed  behind  the  test 
section. 

The  Low  Pressure  System  includes  a  Edwards -Kinney  vacuum  pump,  KD  110,  which  is 
connected  to  the  low  pressure  channel  of  the  shock  tube.  This  tube  can  be 
evaluated  to  less  than  1  mm.  Hg  abs.  The  low  pressure  system  for  the  hypersonic 
nozzle  is  composed  of  three  vacuum  pumps  connected  in  series.  These  include  a 
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rotary  vacuum  pump,  Leybold  S-60,  a  Leybold  Ruvac  25,  of  the  Roots-Blower 
type  and  an  oil  diffusion  pump,  Edward -Speedviac  F-603.  With  this  system 
the  pressure  of  the  nozzle  section  would  be  held  at  approximately  20  y  Hg.  abs. 
The  time  required  to  reach  the  desired  low  pressure  in  the  low  pressure 
channel  was  several  minutes  and  in  the  hypersonic  nozzle,  approximately  15 
minutes . 

The  High  Pressure  Gas  Supply  System  uses  commercial  bottled  high  pressure 
hydrogen  or  the  laboratory  high  pressure- air  supply-  for  air  driven  operation. 
Gas  flow  to  the  shock  tube  is  regulated  by  a  valve  located  on  the  control 
panel,  where  the  compression  chamber  pressure  is  indicated  on  a  pressure  gage. 
This  valve  is  hand  operated  and  is  closed  manually  after  bursting  of  the 
diaphragm. 

4.2.  The  Instrumentation  System 

The  block  diagram  of  the  instrumentation  system  is  presented  in  Fig.  14. 
The  gauges  which  record  the  measured  pressures  prior  to  the  experiment  are 
mounted  on  the  control  panel.  These  measuring  instruments  are:  1.  The  pressure 
in  the  high  pressure-tube-raeasured  by  a  pressure  gauge.  2.  The  low  pressure 
in  the  low  pressure  channel.  3.  The  low  pressure  in  the  supersonic  nozzle. 

4.  Room  Temperature.  The  output  of  the  model  instrumentation  is  recorded  on 
Tetronix  oscilloscopes  models  535,  555  and  565  which  were  operated  in  a  single 
sweep  mode.  The  dr-ta  was  photographed  by  Polaroid  cameras. 
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Pressure  Measurements 

Pressures  were  measured  by  piezoelectric  transducers  manufactured  by  the 
Kistier  Company.  Performance  details  of  the  transducers  are  given  in  the 
following  Table: 


Location  $ 
Type  of 
Measure¬ 
ment 

Model 

No. 

Typical 

Pressure 

Measured 

Range  of 
tranducers 

Elec. 

Output 

Response 

Time 

Natural 

Frequency 

External 
dimensions 
(length  x 
diameter) 

psi 

psi 

pcb/psi 

y  sec 

cps 

inches 

Static 
pressure 
at  the 
end  of 
the  shock 
tube 

601 

300 

10-3000 

■ 

3 

130,000 

0.25  x  0.6 

Total 

pressure 

behind 

the 

shock 

wave  in 

the 

test 

section 

603 

5 

0-3000 

0.37 

1 

400,000 

0.25  x  0.45 

Static 
pressure 
in  the 
test 
section 

606 

0.05 

0-30 

4.4 

3 

130,000 

0.5  x  1.25 

It  is  immediately  obvious  from  the  table  that  except  for  the  transducers  at  the  end  of 
the  shock  tube,  the  others  operated  in  the  lower  portion  of  their  measurement  range. 


this  results  in  a  relatively  low  signal-to-noise  ratio.  In  order  to  improve  the  signal- 
to-noise  ratio  passive  electrical  filters  were  installed.  These  filters  had  two  filter- 
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ing  levels,  one  removes  noise  above  a  frequency  of  1000  Hz  and  the  second 
removes  noise  above  3000  Hz.  The  lag  time  in  the  pressure  measurements  due 
to  these  filters  is  approximately  100  psec.  which  is  tolerable  in  the  shock 
tunnel  measurements.  The  pressure  transducers  were  connected  through  the 
charge  amplifiers  (Kistler  Model  504  amplifiers)  to  the  osciJoscope.  The 
total  pressure  gages  were  mounted  in  a  rake  which  is  shown  in  Fig.  15.  The 
rake  mounting  arrangement  allows  rotation  and  movement  of  the  rake  in  the 
direction  of  the  longitudinal  tunnel  axis.  This  arrangement  enab’es  the 
mapping  of  the  total  pressures  at  various  angular  positions  and  at  various 
axial  positions  of  the  test  section. 

The  static  pressure  was  measured  on  the  test  section  side  wall  through 
specially  made  measuring  ports  in  a  window.  The  pressure  transducer  is  installed 
in  a  plug  which  is  mounted  on  0-rings,  an  arrangements  which  also  provided 
vibration  isolation  of  the  pressure  gages  from  the  tunnel.  Since  all  of  the  plugs 
were  interchangeable,  it  was  possible  to  measure  the  static  pressure  at  various 
positions  along  the  axis  of  the  test  section  and  at  various  heights  in  the  test 
section. 

It  was  found  necessary  to  protect  the  surfaces  of  the  piezoelectric  pressure 
transducers  by  covering  the  face  of  each  gauge  with  a  thin  layer  of  thermal  in¬ 
sulation  material,  such  as  vacuum  grease,  so  that  the  gauge  is  not  exposed  direct¬ 
ly  to  the  hot  ionized  air.  This  protective  coating  eliminated  stray  output  and 
resulted  in  signals  sensitive  to  the  pressure  variations  only.  This  phenomena  is 
mentioned  also  in  Reference  17. 
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Heat  Transfer  Measurements 

Heat  transfer  measurements  were  made  with  the  thin  platinum  film  resistance 
gauges.  These  gauges  are  made  by  sputtering  a  thin  layer  of  platinum  on  pyrex 
glass.  The  output  voltage  variations  as  a  function  of  temperature  change  is 
given  by  the  relation  AV  =  Iq(R  -  Rq)  =  I^aAT,  Response  time  of  these  gauges 
is  smaller  than  a  microsecond.  The  gauges,  their  preparation  and  calibration 
is  described  in  Reference  18.  These  gauges  were  used  in  this  work  for  the 
following  measurements: (1)  measurements  of  the  stagnation  point  heat  transfer 
rates  in  the  hypersonic  flow  in  the  shock  tunnel  test  section  and  (2)  detectors 
of  the  incident  shock  wave  for  the  shock  velocity  measurements  in  the  shock  tube. 

The  gauges  that  were  used  to  measure  heat  transfer  rates  were  sputtered  on 
pyrex  glass  segments  that  were  cut  from  the  hemispherical  head  of  a  15  mm.  dia¬ 
meter  hemispherical  capped  test  tubes.  The  gauges  were  then  inserted  into  steel 
15mm  diameter  hemispherical  recepticles  so  that  smooth  hemisphere  nosed  probes 
were  obtained.  The  sensitive  film  on  the  gauge  was  such  that  the  stagnation  point 
heat  transfer  rate  variation  of  the  portion  of  the  half  sphere  which  it  covered 
was  negligible,  since  the  heat  transfer  rate  near  the  stagnation  points  falls  of 
rather  slowly  (Ref.  19).  The  dimensions  of  the  platinum  film  on  these  gauges  is 
about  5  mm.  long  and  1  mm  wiue,  A  description  of  the  gauge  calibration  and  of  the 
calibration  system  used  in  this  work  may  be  found  in  Ref.  (20).  The  hemispheres 
are  mounted  on  a  rake  which  is  similar  to  the  one  used  for  the  total  pressure 
measurements,  except  for  the  instrumented  probe  section  of  course.  This  rake  is 
shown  in  Fig.  10. 
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The  data  reduction  follows  the  procedures  developed  in  this  laboratory, 
and  others  (References  19  and  20)  to  determine  the  heat  transfer  rates  from 
the  thin  film  gauges  output.  The  electronic  equipment  used  in  this  work 
(Oscilloscopes,  counters,  charge  amplifiers,  pulse  amplifiers  and  electric 
current  supplies)  are  shown  in  Fig.  17. 
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V.  EXPERIMENTAL  RESULTS 
5.1.  Experimental  Results  in  the  Shock  Tube 

The  measurements  in  the  shock  tube  are  compared  with  the  theoretical 
calculations  of  the  flow  in  the  shock  tube  with  reasonably  good  agreement. 

The  incident  shock  wave  attenuation  is  measured  using  four  thin  film  gauges 
positioned  in  the  shock  tube  low  pressure  channel  as  follows:  two  of  the 
gauges  were  located  permanently  at  stations  1  and  2  of  the  shock  tube  (see 
Fig.  18)  and  the  other  two  were  moved  to  the  various  other  measuring  stations. 

It  can  be  seen  in  Fig.  18  that  the  Mach  number  of  the  shock  wave  is  attenuated 
by  1.5%  to  2%  in  the  2  meters  of  the  instrumented  section,  therefore,  for  the 
7  meter  shock  tube  length  the  attenuation  will  be  no  greater  than  5%  to  7%. 

This  shows  that  there  is  only  a  slight  influence  of  the  boundary  layer  growth 
on  the  shock  wave  velocity  in  the  present  shock  tube.  Data  scatter  in  this 
graph  may  be  attributed  to  the  fact  that  of  the  two  counters  that  were  used 
the  Beckman  7360  had  an  accuracy  of  one  usee  and  the  Beckman  6146  had  an 
accuracy  of  0.1  usee.  Since  a  typical  reading  of  the  counters  was  approximate¬ 
ly  150  usee,  then  the  error  caused  by  the  6146  counter  was  negligible  and  that 
caused  by  the  7360  counter  was  about  0.3%. 

Fig.  19  includes  the  compilation  of  results  for  the  relation  between  the 
pressure  ratio  across  the  diaphragm,  p^/p^,  and  incident  shock  wave  Mach  number, 
Mg,  as  measured  between  stations  1  and  2  (see  Fig.  18).  As  seen  in  this  figure, 
the  experimental  data  scatter  is  large.  It  is  possible  to  relate  the  scatter  to 
the  non-uniformity  of  the  diaphragm  material  and  the  non-uniformity  of  the  dia¬ 
phragm  scribbing  process  and  as  a  result  of  this  the  non  uniformity  in  the  dia¬ 
phragm  opening  and  bursting  process.  In  spite  of  the  wide  scatter  it  is  possible 
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to  see  that  the  pressure  ratio  required  in  order  to  obtain  a  given  M  is 
considerably  larger  than  that  predicted  by  the  ideal  gas  theory.  The  pressure 
measurements  at  the  end  of  the  shock  tube  are  presented  in  Figs.  20  and  21, 
and  a  typical  pressure  gauge  output  obtained  on  the  oscilloscope  camera  is 
shown  in  Fig.  22.  It  can  be  seen  that  at  increasing  shock  Mach  numbers  the 

pressure  in  region  2,  P2/P1  increases  above  that  predicted  by  the  ideal  gas 
theory  (Fig.  20).  At  tailoring  conditions,  for  hydrogen  driver  gas,  (i.e.  at 
Ms  s  6.0)  this  deviation  is  still  very  small  (less  than  10%).  The  measure¬ 
ments  of  Pg/Pj  as  a  function  of  are  presented  in  Fig.  21.  It  is  seen  that 

the  experimentally  obtained  values  are  in  good  agreement  (at  least  up  to  the 
tailoring  conditions)  with  the  theoretical  results  for  y  =  1.4,  and  somewhat 
lower  than  those  predicted  for  a  real  gas.  Since  the  measured  p^pj  values 
were  higher  than  the  theoretical  values,  then  the  reflected  shock  wave  must  be 
weaker  than  that  predicted  by  the  theory. 

Fig. 23  shows  the  pressure  of  the  air  at  reservoir  conditions  as  measured  by  the 

pressure  gauge  at  the  end  of  the  shock  tube.  The  lower  trace  corresponds  to  a 
measurement  without  a  nozzle  in  the  high  pressure  region  and  the  upper  one  with 
a  nozzle  there.  From  these  figures  it  is  seen  that  insertion  of  the  nozzle  in 
the  high  pressure  region  results  in  a  lower  pressure  by  preventing  production  of 
the  additional  shock  wave  which  raises  the  pressure.  It  is  seen  that  y>2  am*  P5 
are  equal  in  both  cares;  the  final  pressure,  pg,  of  the  stagnant  air  without  a 
nozzle  in  the  high  pressure  region,  is  approximately  1.6  times  that  which  is 
obtained  with  a  nozzle,  It  should  be  noted  that  the  available  duration  for 
testing  (the  time  during  which  constant  pressure  exists)  is  several  hundred  micro- 


seconds  longer  for  runs  with  the  nozzle  compared  to  those  without  the  nozzle. 

As  the  pressure  gauges  that  were  located  in  the  supersonic  nozzle  were  operat¬ 
ing  at  the  lowest  limit  of  their  range  causing  a  very  low  signal-to-noise 
ratio  and  since  pressure  in  the  supersonic  nozzle  is  a  direct  function  of  pe» 
it  was  preferable  in  most  of  the  present  experiments  to  operate  without  this 
nozzle  in  the  high  pressure  region  in  spite  of  the  fact  that  by  operating  in 
this  manner  the  useful  test  duration  was  reduced  somewhat. 

A  phenomenon  which  was  not  predicted  by  the  theory  and  which  was  present 
in  all  the  experiments  is  seen  well  in  Fig.  23.  This  is  the  momentary  pressure 
drop  immediately  following  the  reflected  shock  wave.  This  drop  lasts  about  500 
nsec  and  afterwards  the  pressure  returns  to  its  previous  value.  In  an  experiment 
in  which  the  end  of  the  shock  tube  was  completely  blocked  the  same  phenomenon  was 
obtained.  Therefore,  the  phenomenon  cannot  be  related  to  the  rarefaction  waves 
produced  with  bursting  of  the  diaphragm  at  the  entrance  to  the  hypersonic  nozzle 
and  to  the  air  flow  into  the  hypersonic  nozzle.  An  explanation  for  this  phenomenon 
which  is  also  mentioned  in  Ref.  (4),  has  not  been  found. 

The  pressure  measurements  enable  as  well  the  evaluation  of  the  lengths  of 
time  of  the  various  processes  in  the  shock  tube  which  can  be  compared  with  the 
theoretically  calculated  processes  as  shown  in  Fig.  4.  The  duration  of  flow 
in  region  2  was  measured  at  the  condition  of  Fig.  22  to  be  43  psec,  which  is  in 
good  agreement  with  the  calculated  value.  The  time  of  arrival  of  the  rarefaction 
wave  for  the  conditions  of  Fig. 23  is  2680  psec  after  the  arrival  of  the  incident 
shock  wave,  This  value  is  larger  than  the  calculated  value  of  2050  psec.  How¬ 
ever,  it  should  be  mentioned  that  the  calculation  of  the  arrival  time  of  the  head 


-  31 


l  i 
£ 


& 

§ 


t 


i 


rarefaction  wave  at  the  right  end  of  the  shock  tube  is  very  inaccurate  since 
the  calculation  neglects  changes  in  thermodynamic  properties,  and  therefore, 
in  the  sonic  velocity  in  the  gas  following  interaction  between  the  shock  and 
rarefaction  waves  travelling  in  the  tube. 

5.2.  Experimental  Results  in  the  Shock  Tunnel 

Calibration  of  the  supersonic  nozzle  is  made  at  two  nozzle  positions: 

★  ★ 
position  A  (area  ratio  A/A  «  225) and  position  B  (area  ratio  A/A  =  500) 

(see  Fig.  25) .  At  each  nozzle  setting  the  total  pressures  in  the  test  section 

is  measured  with  the  pressure  rake  and  the  stagnation  point  heat  transfer  rate 

with  the  thin  film  resistance  gages  rake.  The  measurements  are  done  at  two 

stations  on  the  nozzle  longitudinal  axis,  stations  1  and  2,  as  shown  in  Fig. 

25.  At  each  of  these  stations  horizontal  and  vertical  distributions  are  measured. 

Static  pressure  gauges  are  positioned  at  various  points  on  the  sidewall  window 

of  the  test  section,  and  both  axial  and  vertical  static  pressure  distributions 

in  the  test  section  were  measured. 


The  pressure  gauge  output  is  determined  directly  from  the  photographs  of 
the  pressure  traces  by  measurement  of  the  difference  between  the  pressure  during 
the  experiment  and  the  initial  pressure.  From  knowledge  of  the  nozzle  initial 
pressure  it  is  possible  to  determine  the  value  of  the  absolute  pressure  during 
the  experiment. 


The  heat  transfer  rate  is  determined  using  the  following  relation: 

.  r  ,  ,0.5  AT  ir  (pCpK)  *  ,1  ,  AV 

q  =  ft!  (pL  h),  ]  — -  *  rr  - c -  (nT1-)  

H  p  bJ  J?  2  a  Vo 
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The  parameter  - =•  is  the  calibration  coefficient  of  the  gauge  which 

<’Cpk)b 

is  obtained  during  its  calibration,  Iq  and  Rq  are  measured,  and  then 

q  is  determined.  In  this  work  a  computer  program  was  used  to  find  the 

best  fit  parabola  for  the  measured  values  of  AV  vs.  t,  thus  determining 
AV 
t 

Fig.  26a  shows  the  trace  of  the  total  pressure  behind  the  shock  wave 
for  a  pressure  ratio  p^/p^  =  610.  This  pressure  ratio  gives  an  indication 
of  the  low  pressure  in  the  tunnel  prior  to  the  experiment  since  the  pressure 
Pj  was  approximately  35  mm  Hg  for  all  of  the  experiments,  The  expansion  fan 
produced  after  bursting  of  the  cellophane  diaphragm  is  seen  in  this  figure. 

Fig.  26b  also  shows  the  trace  of  the  total  pressure  behind  the  shock  wave,  how¬ 
ever,  in  this  case  the  pressure  ratio  was  P1/pi  *  185.  In  this  case  it  is  seen 
that  the  expansion  fan  at  the  beginning  of  the  flow  passes  through  a  shock  wave 
which  cause  a  momentary  pressure  rise.  Fig.  26c  shows  the  static  pressure  trace. 
Typical  output  of  the  stagnation  point  heat  transfer  rate  are  shown  in  Fig.  27. 

The  influence  of  the  starting  phenomenon  on  the  test  duration  is  demonstrated  in 
Fig.  28.  Since  the  static  pressure  for  the  nozzle  in  position  B  is  considerably 
lower  than  that  of  position  A  the  duration  of  the  starting  process  in  this  nozzle 
at  position  A  is  much  shorter  and  steady  conditions  are  obtained  from  t  *  1.8  to 
2.9  msec,  while  with  the  nozzle  at  position  B  the  steady  conditions  are  obtained 
only  after  2.4msec  t.o  3.0msec.  These  traces  are  typical  of  all  the  experiments  that 
were  performed  in  the  hypersonic  nozzle. 

Figs.  29  to  32  summarize  the  total  pressure  distribution  measured  in  the 
test  section.  It  is  seen  that  in  both  nozzle  positions  the  horizontal  pressure 
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distribution  is  uniform  and  the  boundary  layer  is  relatively  thin  (Figs. 

29  and  31).  When  the  nozzle  is  in  position  A,  two  pressure  levels  appear 
in  the  vertical  distribution  (Fig.  30),  in  the  lower  half  of  the  cross- 
section  the  total  pressure  behind  the  shock  wave  is  higher  than  that  in  the 
upper  half.  The  boundary  layer  in  this  case  is  also  relatively  thin.  The 
vertical  mapping  with  the  nozzle  in  position  B(Fig.  32)  shows  a  uniform 
pressure  distribution  in  the  center  of  the  cross-section,  however,  the 
pressure  drop  away  from  the  center  core  indicates  a  relatively  thick  boundary 
layer. 

Vertical  mappings  of  the  static  pressure  in  the  test  section  (Figs.  33 
and  34)  verify  the  phenomenon  observed  with  the  total  pressure  gauges,  The 
static  pressure  mapping  of  position  A,  Fig.  33,  also  shows  two  pressure 
levels,  the  higher  level  being  located  in  the  lower  half  of  the  section.  Fig. 
34  shows  a  uniform  vertical  static  pressure  distribution  when  the  nozzle  is 
in  position  B.  The  stagnation  point  heat  transfer  distribution  in  the  test 
section,  shown  in  Figures  35  to  38,  is  very  similar  to  the  pressure 
distributions.  The  horizontal  heat  transfer  distributions  are  uniform  for 
both  nozzle  position,  as  seen  in  Figs.  35  and  37.  In  both  of  these  nozzle 
positions  there  is  an  indication  of  a  thin  boundary  layer,  In  the  vertical 
heat  transfer  rate  distributions  for  the  nozzle  position  A,  two  heat  transfer 
rate  levels  are  found  where  the  higher  level  is  found  in  the  lower  part  of  the 
section  (  although  the  phenomenon  is  smaller  than  in  the  pressure  distribution 
case) .  The  uniform  core  at  the  center  of  the  test  section  is  narrower  in  the 
heat  transfer  measurements  than  that  obtained  from  the  total  pressure  measure¬ 


ments  . 
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The  static  and  total  pressure  measurements  enable  the  determination  of  the 
test  section  Mach  number  distribution.  In  this  case  an  equilibrium  expansion  in 
the  nozzle  was  assumeu  for  the  Mach  number  evaluation. 

For  nozzle  position  B  there  is  very  good  agreement  between  the  Mach  number 
that  was  obtained  from  the  total  pressure  measurements  behind  the  shock  wave 
and  that  obtained  from  the  static  pressure  measurements, (M  =  6.5?  to  7.0  was 
obtained) . 

The  nozzle  at  position  A  has  a  somewhat  non  uniform  flow  due  probably 
to  various  disturbances  which  are  due  to  the  asymetric  wedge  type  nozzle  shape. 
Therefore  for  the  nozzle  in  position  A  the  Mach-  numbers  obtained  from  the  total 
pressure  measurements  are  different  from  those  obtained  from  the  static  pressure 
measurements.  These  differences  indicate  a  Mach  number  of  about  5.1  in  the 
lower  part  of  the  nozzle  and  5.7  in  the  upper  part. 

Good  agreement  is  obtained  when  the  heat  transfer  rates  that  were  obtained 

in  the  experiments  (Figs.  35  to  38)  are  compared,  with  the  theoretical  results 

(Figs.  10  and  10A)  when  the  flew  Mach  number  is  the  one  that  is  determined  from 

the  pressure  measurements  and  the  velocity  gradient  at  the  stagnation  point 
dU 

(-jp) s  is  calculated  using  the  potential  flow  distribution  arour.d  a  sphere  (see 
Appendix  A) .  When  this  velocity  gradient  term  is  calculated  according  to  the 
Newtonian  approximation  (also  found  in  Appendix  A),  the  experimental  results 
are  40%  lower  than  the  theoretical  results  at  the  assumed  flow  Mach  number. 
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VI •  SUMMARY  AND  CONCLUSIONS 

The  present  shock  tunnel  was  added  to  an  existing  shock  tube  in  which 
the  ratio  between  the  lengths  of  the  driver  gas  section  and  the  driven  gas 
sectaai  is  suitable  for  work  in  the  shock  tube  itself  since  the  driven  gas 
section  is  long  md  gives  a  relatively  long  time  between  the  incident  shock 
wave  and  the  contact  surface  (Region  2).  However,  this  length  ratio  is  not 
the  optimal  ratio  for  maximum  flow  duration  in  the  hypersonic  nozzle  for  the 
hydrogen/air  tailoring  conditions.  It  is  shown  that  by  changing  this  length 
ratio  it  will  be  possible  to  increase  the  flow  duration  in  the  hypersonic 
nozzle  by  a  factor  of  two.  It  is  also  shown  that  high  initial  pressures  are 
very  advantageous  since  they  enable  reasonable  pressures  in  the  hypersonic 
nozzle  test  section.  The  extremely  simple  nozzle  design  used  in  the  present 
shock  tunnel  enabled  easy  regulation  of  the  flow  Mach  number  in  the  test 
section  and  surprisingly  good  flow  properties. 

This  work  indicates  that  it  is  possible  to  use  the  present  shock  tunnel 
for  aerodynamic  measurements  at  reasonably  high  stagnation  enthalpies  (2000 
Btu/lb.  for  hydrogen/air  at  tailoring  conditions)  and  high  Mach  numbers. 

Test  time  of  steady  flow  duration  in  the  hypersonic  nozzle  is  at  least  SOOusec, 
for  the  Mach  7  nozzle  and  about  lysec.  for  the  Mach  5.5  nozzle.  A  more  optimal 
design  will  result  in  about  2  msec,  of  steady  flow  in  the  hypersonic  nozzle. 

The  uniform  flow  core  in  the  nozzle  is  about  8  cm  high  and  about  16  cm  wide, 
so  that  relatively  large  models  can  be  used  in  this  test  section  for  hyper¬ 
sonic  testing. 
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The  average  properties  in  the  test  section  for  the  two  nozzle  positions 
determined  in  the  present  investigation  are: 


Nozzle 

Position 

Nominal 

Mach 

Number 

P 

ata 

(mm  Mg) 

Po2 

ata 

B 

U 

m/sec 

Re/m 

1/m 

A 

5.5 

6.8  x  10‘3 
(5) 

0.33 

550 

2850 

3.6  x  105 

B 

? 

2.5  x  10"3 
(2) 

0.16 

420 

2900 

2.1  x  105 
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APPENDIX  A 

CALCULATION  OF  THE  VELOCITY  GRADIENT  AT  THE  STAGNATION  POINT  OF  A  SPHERE 

IN  HYPERSONIC  FLOW 

Since  the  flow  behind  the  strong  detached  shock  wave  is  low  subsonic 
flow,  it  is  possible  to  assume  that  the  density  at  the  stagnation  region 
is  constant  and  equal  the  density  at  the  stagnation  point.  In  this  case 
the  velocity  distribution  in  the  stagnation  region  can  be  calculated  from 
the  subsonic  incompressible  potential  solution  of  the  flow  over  a  sphere. 

The  potential  equation  is  then 

Ur^ 

U  2  .  2  o  .2 

<J>  a  -  y  r  sin  e  +  2-f~  sin  0 


where: 

Ur  *  radial  velocity 

U.  »  velocity  in  direction  0 
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On  the  sphere  surface,  r  =  rQ,  Ur  *  0  and  *  3/2  U  sine 
Since 

d  =  r  d0 
s  o 


and 


dU 


dU 


ds  "  de 

at  the  stagnation  point 


de 

ds 


2r 


U  cos  0 


A 

Ws 


From  this  relation: 


U 


U  = 


-  P) 


0.5 

■1 


where: 

p  -  static  pressure  behind  a  shock  wave 
pg  -  stagnation  pressure  behind  a  shock  wave. 


So  that: 


P  “  P  0.5 

©s  ■  5ri2i— 1 

0  s 


Pc-  P  °'5 

7-  t-f — 1 

2  r  s 
0 


The  velocity  gradient  can  be  also  estimated  assuming  the  Newtonian  pressure 


distribution  as  follows: 


C 

Pmax 


P  -  P„ 


1/2  p  tr 
00  00 

PSJ  p» 

1/2  77 

00  CO 


C  (1  -  sin20) 
Pmax 


where: 

p  -  static  pressure  on  the  sphere  surface,  static 
pressure  behind  a  shock  wave. 

ps  -  stagnation  pressure  behind  a  shock  wave 

p^  -  freestream  pressure  in  front  of  a  shock  wave 

-  freestream  velocity  in  front  of  a  shock  wave 


and  the  velocity  gradient  at  the  stagnation  point  is  in  this  case: 
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FIG.  NO.  22:  TRACES  OF  PRESSURE  MEASUREMENT  AT  THE  END 
OF  THE  SHOCK  TUBE. 
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FIG.  NO.  23:  TRACES  OF  PRESSURE  MEASUREMENT  AT  THE  END 
OF  THE  SHOCK  TUBE. 
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FIG.  NO.  24:  RATIO  OF  EQUILIBRIUM  PRESSURE  AT  THE  END  OF  THE 
SHOCK  TUBE  TO  THE  INITIAL  PRESSURE  IN  THE  SHOOK 
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FIG. NO.  26:  PRESSURE  GAUGE  TRACES  IN  TEST  SECTION. 
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FIG.  NO.  27:  HEAT  TRANSFER  RATE  MEASUREMENTS  TRACES  IN 
TEST  SECTION. 
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OKTaL  DISTRIBUTION  OF  STAGNATION  POINT  HEAT  TRANSFER  RATES  NOZZLE  IN 
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